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We study dynamics of an antiferromagnetic soliton under a temperature gradient. To this end, 
we start by phenomenologically constructing the stochastic Landau-Lifshitz-Gilbert equation for an 
antiferromagnet with the aid of the fluctuation-dissipation theorem. We then derive the Langevin 
equation for the soliton’s center of mass by the collective coordinate approach. An antiferromagentic 
soliton behaves as a classical massive particle immersed in a viscous medium. By considering a 
thermodynamic ensemble of solitons, we obtain the Fokker-Planck equation, from which we extract 
the average drift velocity of a soliton. The diffusion coefficient is inversely proportional to a small 
damping constant a, which can yield a drift velocity of tens of m/s under a temperature gradient 
of 1 K/mm for a domain wall in an easy-axis antiferromagnetic wire with a ~ 10~^. 
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Introduction .—Ordered magnetic materials exhibit 
solitons and defects that are stable for topological rea¬ 
sons [T]. Well-known examples are a domain wall (DW) 
in an easy-axis magnet or a vortex in a thin film. Their 
dynamics have been extensively studied because of fun¬ 
damental interest as well as practical considerations such 
as the racetrack memory [2] . A ferromagnetic (FM) soli¬ 
ton can be driven by various means, e.g., an external 
magnetic field [3] or a spin-polarized electric current [4]. 
Recently, the motion of an FM soliton under a temper¬ 
ature gradient has attracted a lot of attention owing to 
its applicability in an FM insulator [5]-[8]. A temperature 
gradient of 20 K/mm has been demonstrated to drive a 
DW at a velocity of 200 /im/s in an yttrium iron garnet 
film [9]. 

An antiferromagnet (AFM) is of a great current inter¬ 
est in the field of spintronics [TOHE] due to a few advan¬ 
tages over an FM. First, the characteristic frequency of 
an AFM is several orders higher than that of a typical 
FM, e.g., a timescale of optical magnetization switching 
is an order of ps for AFM NiO [13] and ns for FM Cr02 
m, which can be exploited to develop faster spintronic 
devices. Second, absence of net magnetization renders 
the interaction between AFM particles weak, and, thus, 
leads us to prospect for high-density AFM-based devices. 
Dynamics of an AFM soliton can be induced by an elec¬ 
tric current or a spin wave UMI]. 

A particle immersed in a viscous medium exhibits a 
Brownian motion due to a random force that is required 
to exist to comply with the fluctuation-dissipation theo¬ 
rem (FDT) [T9] . An externally applied temperature 
gradient can also be a driving force, engendering a phe¬ 
nomenon known as thermophoresis [20] . Dynamics of an 
FM and an AFM includes spin damping, and, thus, in¬ 
volves thermal fluctuations at a finite temperature [21]. 
The corresponding thermal stochastic field influences dy¬ 
namics of a magnetic soliton [Hl[22l|23], e.g., by assisting 
a current-induced motion of an FM DW [24] . 


Brownian motion of a DW 



FIG. 1. (Golor online) A thermal stochastic force caused by 
a temperature gradient pushes an antiferromagnetic domain 
wall to a colder region. The diffusion coefficient of the domain 
wall is inversely proportional to a small damping constant, 
which may give rise to a sizable drift velocity. 

In this Rapid Communication, we study the Brown¬ 
ian motion of a soliton in an AFM under a tempera¬ 
ture gradient. We derive the stochastic Landau-Lifshitz- 
Gilbert (LLG) equation for an AFM with the aid of the 
FDT, which relates the fluctuation of the staggered and 
net magnetization to spin damping. We then derive the 
Langevin equation for the soliton’s center of mass by em¬ 
ploying the collective coordinate approach [Hes]. We 
develop the Hamiltonian mechanics for collective coordi¬ 
nates and conjugate momenta of a soliton, which sheds 
light on stochastic dynamics of an AFM soliton; it can be 
considered as a classical massive particle moving in a vis¬ 
cous medium. By considering a thermodynamic ensemble 
of solitons, we obtain the Fokker-Planck equation, from 
which we extract the average drift velocity. As a case 
study, we compute the drift velocity of a DW in a quasi 
one-dimensional easy-axis AFM. 

Thermophoresis of a Brownian particle is a multi¬ 
faceted phenomenon, which involves several competing 
mechanisms. As a result, a motion of a particle depends 
on properties of its environment such as a medium or 
a temperature T [26]. For example, particles in pro¬ 
tein (e.g., lysozyme) solutions move to a colder region 
for T > 294 K and otherwise to a hotter region ISniET]. 
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Thermophoresis of an AFM soliton would be at least as 
complex as that of a Brownian particle. We focus on one 
aspect of it in this Rapid Communication; the effect of 
thermal stochastic force on dynamics of the soliton. We 
discuss two other possible mechanisms, the effects of a 
thermal magnon current and an entropic force [5], later 
in the Rapid Communication. 

Main results. —Before pursuing details of derivations, 
we first outline our three main results. Let us consider 
a bipartite AFM with two sublattices that can be trans¬ 
formed into each other by a symmetry transformation 
of the crystal. Its low-energy dynamics can be devel¬ 
oped in terms of two fields: the unit staggered spin 
field n = (mi — m 2)/2 and the small net spin field 
m = (mi + m 2)/2 perpendicular to n. Here, mi and 
m 2 are unit vectors along the directions of spin angular 
momentum in the sublattices. 

Starting from the standard Lagrangian description of 
the antiferromagnetic dynamics [28], we will show below 
that the appropriate theory of dissipative dynamics of 
antiferromagnets at a finite temperature is captured by 
the stochastic LLG equation 

s{h -h Pn X ih) =n x (h + h^^), (la) 

s{ih -h ydm X ih -h an x h) =n x (g + 

+ mx (h + h^^), (lb) 

in conjunction with the correlators of the thermal 
stochastic fields g^^ and 

t')) = 2kBTas5ij5{r - r')5{t - t'), (2a) 

,t')) = 2kBTl3s5ij5{r — r')5{t — t'), (2b) 

which are independent of each other [29] . This is our first 
main result. Here, a and P are the damping constants 
associated with h and ih, g = —5U/5n and h = —5U/5m 
are the effective fields conjugate to n and m, l/[n, m] = 
l/[n] + j dV\mp /2x is the potential energy (y represents 
the magnetic susceptibility), and s = hS/V is the spin 
angular momentum density (V is the volume per spin) 
per each sublattice. The potential energy l/[n(r,t)] is 
a general functional of n, which includes the exchange 
energy j dVAijdin • djn at a minimum [28] . 

Slow dynamics of stable magnetic solitons can often 
be expressed in terms of a few collective coordinates 
parametrizing slow modes of the system. The center of 
mass R represents the proper slow modes of a rigid soli¬ 
ton when the translational symmetry is weakly broken. 
Translation of the stochastic LLG equation 0 into the 
language of the collective coordinates results in our sec¬ 
ond main result, a Langevin equation for the soliton’s 
center of mass R: 

MR + TK = -dU/dR + F^'^, (3) 

which adds the stochastic force F**’ to Eq. (5) of Tveten 
et al. m- The mass and dissipation tensors are symmet¬ 
ric and proportional to each other: Mij = p f dV{din • 


djn) and Tij = Mijjr, where r = p/as is the relaxation 
time, p = is the inertia of the staggered spin field 
n. The correlator of the stochastic field obeys the 
Einstein relation 

= 2kBTTij5(t - tj. (4) 

A temperature gradient causes a Brownian motion of 
an AFM soliton toward a colder region. In the absence 
of a deterministic force, the average drift velocity is pro¬ 
portional to a temperature gradient V ex ksVT in the 
linear response regime. The form of the proportional¬ 
ity constant can be obtained by a dimensional analy¬ 
sis. Let us suppose that the mass and dissipation ten¬ 
sors are isotropic. The Langevin equation 0 is, then, 
characterized by three scalar quantities: the mass M, 
the viscous coefficient F, and the temperature T, which 
define the unique set of natural scales of time r = M/F, 
length I = y/ksTM/T ^ and energy e = ksT. Using 
these scales to match the dimension of a velocity yields 
V = —cp{kB'VT), where p = F”^ is the mobility of an 
AFM soliton and c is a numerical constant. The explicit 
solution of the Fokker-Planck equation, indeed, shows 
c = 1. This simple case illustrates our last main result; 
a drift velocity of an AFM soliton under a temperature 
gradient in the presence of a deterministic force F is given 
by 

V = pF - p{kBVT). ( 5 ) 

For a DW in an easy-axis one-dimensional AFM, the 
mobility is p = X/2asa, where A is the width of the 
wall and a is the cross-sectional area of the AFM. For 
a numerical estimate, let us take an angular momen¬ 
tum density s = 2/inm“^, a width A = 100 nm, and a 
damping constant a = 10“^ following the previous stud¬ 
ies [niisn]. For these parameters, the AFM DW moves 
at a velocity V = 32 m/s for the temperature gradient of 
VT = 1 K/mm. 

Stochastic LLG equation. —Long-wave dynamics of an 
AFM on a bipartite lattice at zero temperature can de¬ 
scribed by the Lagrangian [28] 

L = s J dVm • (n X h) — l/[n, m]. (6) 

We use the potential energy l/[n, m] = f dU|mp/2y + 
U[n] throughout the Rapid Gommunication, which re¬ 
spects the sublattice exchange symmetry (n — n, m ^ 
m). Minimization of the action subject to nonlinear con¬ 
straints |n| = 1 and n- m = 0 yields the equations of mo¬ 
tion for the fields n and m. Damping terms that break 
the time reversal symmetry can be added to the equa¬ 
tions of motion to the lowest order, which are first order 
in time derivative and zeroth order in spatial derivative. 
The resultant phenomenological LLG equations are given 
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by 


s{h + X ih) = n X h, (7a) 

s(m + f3m X ih + an xh) = nxgH-mxh (7b) 

[UlsollSI]. The damping terms can be derived from the 
Rayleigh dissipation function 

R = y* dV{as\nf l3s\mf)/2, (8) 

which is related to the energy dissipation rate hy —U = 
2R. The microscopic origin of damping terms does not 
concern us here but it could be, e.g., caused by thermal 
phonons that deform the exchange and anisotropy inter¬ 
action. 

At a finite temperature, thermal agitation causes fluc¬ 
tuations of the spin fields n and m. These thermal fluc¬ 
tuations can be considered to be caused by the stochas¬ 
tic fields and with zero mean, which are con¬ 
jugate to n and m, respectively; their noise correlators 
are then related to the damping coefficients by the FDT. 
The standard procedure to construct the noise sources 
yields the stochastic LLG equation 0- The correlators 
of the stochastic fields are obtained in the following way 
[T8l[32]. Casting the linearized LLG equation 0 into the 
form {h, g} = 7(8) {h, ih} provides the kinetic coefficients 
7. Symmetrizing the kinetic coefficients 7 produces the 
correlators Q of the stochastic fields consistent with the 
FDT. 

Langevin equation. —For slow dynamics of an AFM, 
the energy is mostly dissipated through the temporal 
variation of the staggered spin field n due to |mp 
(ar)^|hp ^ |hp (from Eq. Q), which allows us to set 
/3 = 0 to study long-term dynamics of the magnetic soli- 
ton m- At this point, we switch to the Hamiltonian for¬ 
malism of an AFM [33] , which sheds light on the stochas¬ 
tic dynamics of a soliton. The canonical momentum field 
TT conjugate to the staggered spin field n is 

TT = 6 L/ 6 h = 5m X n. ( 9 ) 

The stochastic LLG equations 0 can be interpreted as 
Hamilton’s equations, 

h = 6 HI 67 V = TV I p, TT = — 6 H/ 6 n — 6 R/ 6 h-\-g^^^ (10) 
with the Hamiltonian 

H = J dVTT-h~L = j dV^-^+U[n]. ( 11 ) 

Long-time dynamics of magnetic texture can often be 
captured by focusing on a small subset of slow modes, 
which are parametrized by the collective coordinates 
q = {^1: ^2, • • • }• A classical example is a DW in a one¬ 
dimensional easy-axis magnet described by the position 
of the wall X and the azimuthal angle ^ jS] |33]. An¬ 
other example is a skyrmion in an easy-axis AFM film. 


which is described by the position R = {X^Y) [34l [35] . 
Translation from the field language into that of collective 
coordinates can be done as follows. If the staggered spin 
field n is encoded by coordinates q as n(r, t) = n[r; q(t)], 
time dependence of n reflects evolution of the coordi¬ 
nates: h = qidnjdqi. With the canonical momenta p 
defined by 


= [ 


, dn 


Pi = ^ 

Hamilton’s equations (10) translate into 

Mci = p, p + rq = F + F*'', 


( 12 ) 


(13) 


where F = —dU/dq is the deterministic force and = 
f dVdq. n ■ g**’ is the stochastic force. Hamilton’s equa- 
tions (H^ can be derived from the Hamiltonian in the 
collective coordinates and conjugate momenta. 


H = pTM-ip/2 + C/(q), (14) 


with the Poisson brackets {qi^Pj} = dij, {qi^qj} = 
{Pi^Pj} = 0- AFM soliton, thus, behaves as a classi¬ 
cal particle moving in a viscous medium. 

We focus on a translational motion of a rigid AFM 
soliton by choosing its center of mass as the collective 
coordinates q = R; n(r, t) = n(r — R(t)). Eliminat¬ 
ing momenta from Hamilton’s equations (13) yields the 
Langevin equation for the soliton’s center of mass: 


tR T R — pF T ?7, 

where r] = is the stochastic velocity. Here the mo¬ 
bility tensor of the soliton /r = F”^ relates a deterministic 
force to a drift velocity (R) = pF at a constant temper¬ 
ature [36|. The mobility is inversely proportional to a 
damping constant, which can be a small number for an 
AFM, e.g., a ^ 10“^ for NiO [37|. The correlator ([^ of 
thermal stochastic fields is translated into the correlator 
of the stochastic velocity. 


{r]i{t)rjj{t')) = 2 kBTpij 6 {t - t') = 2 Dij 6 {t - t'). (16) 


From Eq. (16), we see that diffusion coefficient and the 


mobility of the soliton respect the Einstein-Smoluchowski 
relation: D = pksT^ which is expected on general 
grounds. It can also be explicitly verified as follows. 
A system of an ensemble of magnetic solitons at ther¬ 
mal equilibrium is described by the partition function 
Z = j Iii[dpidxi/2'Kh\ exp(—i7//cBT), which provides the 
autocorrelation of the velocity, {xiXj)l 2 = M-^^/cbT/ 2 
(the equipartition theorem). In the absence of an ex¬ 
ternal force, multiplying rxi Y Xi = pi ( [T^ by Xj and 
symmetrizing it with respect to indices i and j give the 
equation, Td‘^{xiXj)ldt‘^ + d{xiXj)/dt = 2r(iiXj), where 
the first term can be neglected for long-term dynamics 
This equation in conjunction with the autocorre¬ 
lation of the velocity allows us to obtain the diffusion co¬ 


efficient Dij in Eq. (16), (xiXj) = 2kBTTM^-^t = 2Dijt, 
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without prior knowledge about the correlator (§ of the 
stochastic fields. 

Average dynamies .—An AFM soliton exhibits Brown¬ 
ian motion at a finite temperature. The following Fokker- 
Planck equation for an ensemble of solitons in an inho¬ 
mogeneous medium describes the evolution of the density 
p(R, t) at time t ^ r: 

^ + V-j = 0, with j = iJ,Fp-DVp-DT{kB'VT), (17) 

where Dt = jap is the thermophoretic mobility (also 
known as the thermal diffusion coefficient) [20l [38]. A 
steady-state current density j = /iFpo — VT) with 

a constant soliton density p(r, t) = po solves the Fokker- 
Planck equation from which the average drift veloc¬ 
ity of a soliton can be extracted [39] : 

V = pF - p{kBVT). (18) 


Let us take an example of a DW in a quasi one¬ 
dimensional easy-axis AFM with the energy U[n] = 
J dV{A\dxn\‘^ — Kn‘l)l2. A DW in the equilibrium 
is = (sin^ cos 4>, sin^ sin <F, cos ^) with cos^ = 

tanh[(x — X)/A], where A = yjAjK is the width of 
the wall. The position X and the azimuthal angle 4> 
parametrize zero-energy modes of the DW, which are en¬ 
gendered by the translational and spin-rotational symme¬ 
try of the system. Their dynamics are decoupled, Tx^ = 
0, which allows us to study the dynamics of X separately 
from 4>. The mobility of the DW is p = A/2(ascr, where a 
is the cross-sectional area of the AFM. The average drift 
velocity (18) is given by 


y = 


1 k^AXT 
2 a S(j 


(19) 


Diseussion —The deterministic force F on an AFM 
soliton can be extended to include the effect of an elec¬ 
tric current, an external field, and a spin wave [IMZ]. It 
depends on details of interaction between the soliton and 
the external degrees of freedom, whose thorough under¬ 
standing would be necessary for a quantitative theory for 
the deterministic drift velocity pF. The Brownian drift 
velocity V ( [T8| ) is, however, determined by local property 
of the soliton. We have focused on the thermal stochastic 
force as a trigger of thermophoresis of an AFM soliton in 
this Rapid Communication. There are two other possi¬ 
ble ingredients of thermophoresis of a magnetic soliton. 
One is a thermal magnon current, scattering with which 
could exert a force on a soliton [40]. The other is an 
entropic force, which originates from thermal softening 
of the order-parameter stiffness [5]. Effects of these two 
mechanisms have not been studied for an AFM soliton; 
full understanding of its thermophoresis is an open prob¬ 
lem. 

In order to compare different mechanisms of thermally- 
driven magnetic soliton motion, let us address a closely 


related problem of thermophoresis of a DW in a quasi 
one-dimensional FM wire with an easy-xz-plane easy- 
2 )-axis |3], which has attracted a considerable scrutiny 
recently. To that end, we have adapted the approach de¬ 
veloped in this Rapid Communication to the FM case, 
which leads to the conclusion that a DW drifts to a 
colder region by a Brownian stochastic force at the ve¬ 
locity given by the same expression for an AFM DW, 
= —kBX^T/2as(j [41]. A thermal magnon cur¬ 
rent pushes a DW to a hotter region at the velocity 
= kBXT/67r‘^sXm, where A^ = ^/HA/sT is the 
thermal-magnon wavelength [7]. According to Schlick- 
eiser et al [5] , an entropic force drives a DW to a hotter 
region at the velocity = /c^VT/dsa, where a is the 
lattice constant. The Brownian stochastic force, there¬ 
fore, dominates the other forces for a thin wire, a ^ Xa/a 
(supposing rigid motion) [42] . 

Within the framework of the LLG equations that are 
first order in time derivative, the thermal noise is white 
as long as slow dynamics of a soliton is concerned, 
i.e., the highest characteristic frequency of the natural 
modes parametrized by the collective coordinates is much 
smaller than the temperature scale, liw kBT. The 
thermal noise could be colored in general [26], e.g., for 
fast excitations of magnetic systems, which may be ex¬ 
amined in the future. In addition, local energy dissi¬ 
pation © allowed us to invoke the standard FDT at 
the equilibrium to derive the stochastic fields. It would 
be worth pursuing to understand dissipative dynamics 
of general magnetic systems, e.g., with nonlocal energy 
dissipation with the aid of generalized FDTs at the out- 
of-equilibrium [43] . 


We have studied dynamics of an AFM soliton in the 
Hamiltonian formalism. Hamiltonian’s equations (13) for 
the collective coordinates and conjugate momenta can be 
derived from the Hamiltonian (14) with the conventional 
Poisson bracket structure. By replacing Poisson brackets 
with commutators, the coordinates and conjugate mo¬ 
menta can be promoted to quantum operators. This may 
provide a one route to study the effect of quantum fluc¬ 
tuations on dynamics of an AFM soliton [44]. 


After the completion of this work, we became aware 
of two recent reports. One is on thermophoresis of an 
AFM skyrmion [35], whose numerical simulations sup¬ 
port our result on diffusion coefficient. The other is 
on thermophoresis of an FM DW by a thermodynamic 
magnon recoil [45] . 
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